Plasma production and thermalisation in a strong field 
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Aspects of the formation and equilibration of a quark-gluon plasma are explored using a 
quantum kinetic equation, which involves a non-Markovian, Abelian source term for quark and 
antiquark production and, for the collision term, a relaxation time approximation that defines 
a time-dependent quasi-equilibrium temperature and collective velocity. The strong Abelian 
field is determined via the simultaneous solution of Maxwell's equation. A particular feature of 
this approach is the appearance of plasma oscillations in all thermodynamic observables. Their 
presence can lead to a sharp increase in the time-integrated dilepton yield. 
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I. INTRODUCTION 

The Relativistic Heavy-Ion Collider (RHIC) at Brook- 
haven National Laboratory and the Large Hadron Col- 
hder (LHC) at CERN are designed with the goal of 
producing an equilibrated phase of deconfined partonic 
matter: the quark gluon plasma (QGP). Lattice-QCD 
simulations 111 and well-constrained phenomenological 
models ||^ predict a second order phase transition at 
Tc — 150 MeV; i.e., at energy densities ^ IGeV/fm^, in 
the equilibrium two light-flavour theory. For more than 
two flavours the characteristics of the transition are not 
as clear, but there is a transition. A number of phenom- 
ena have been proposed as signals for the existence of an 
equilibrated QGP |^ but the nonequilibrium stages in 
the plasma's development are poorly understood and it 
is a contemporary challenge to develop a description of 
the spacetime evolution of an ultra-relativistic heavy-ion 
collision (URHIC): from particle production in the col- 
lision, through equilibration and plasma formation, and 
on to hadronisation. 

Two methods are commonly used to describe the pro- 
duction of partons in a collision: a perturbative pre- 
formed parton picture Q and a nonperturbative flux- 
tube-based picture j^. They are complementary, and 
Monte-Carlo event generators 0, || and hydrodynam- 
ical models |^ have been developed to facilitate the anal- 
ysis of data using either production model. 

Herein we focus on the pre-equilibrium particle- 
production stage in the evolution of a QGP and choose to 
employ a flux-tube model. In this model the two collid- 
ing nuclei are imagined to pass through one another and 
stretch a high energy-density flux-tube between them- 



selves as they separate. This flux tube, which describes 
the highly excited QCD vacuum, decays via a nonpertur- 
bative particle-antiparticle production process analogous 
to the Schwinger mechanism. 

Particle production by the flux tube is described by a 
source term and in quantum flcld theory that source term 
is non-Markovian 0, |ll|, |l|, |l3[ |l|, |^ ; i.e., essentially 
nonlocal in time. This feature can be important when 
the flelds are strong. (In the weak-field limit a time-local 
Schwinger-like source term is recovered.) Such a situa- 
tion is plausible at RHIC and especially at LHC |l^ 
where the anticipated initial energy densities are, respec- 
tively, e ~ 10-100 GcV/fm3 and £ > 1 TeV/fm^. 

Another feature that is characteristic of the flux-tube 
production mechanism is the back-reaction phenomenon. 
This phrase simply describes the fact that once the par- 
ticles are produced they are accelerated and thereby gen- 
erate a field that interferes with the coUisional field that 
produced them. Plasma oscillations are then almost 
inevitable p^ , although they can rapidly be damped 
if the thermalising collisions between particles are fre- 
quent 0, |l| |0[ |l], H, g. 

An observable that may preserve information about 
non-Markovian effects and plasma oscillations in the pre- 
equilibrium particle-production stage of an URHIC is the 
thermal dilepton spectrum because leptons do not par- 
ticipate in the strong interactions that equilibrate the 
QGP iQ. To explore that possibility, we employ a quan- 
tum Vlasov equation with a non-Markovian source term 
to calculate the single particle distribution function that 
characterises particle production by the flux tube, /(p, t), 
and use that to calculate the dilepton spectrum and its 
evolution from impact to equilibration. 
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Of course, equilibration can only be effected by dissipa- 
tive processes, such as collisions, and herein we describe 
those effects via a relaxation time approximation (RTA) . 
This is a coarse representation of the interactions be- 
tween the partons produced in the collision but, even so, 
ensuring thermodynamic consistency is nontrivial. We 
introduce and describe one practical scheme for achiev- 
ing that goal. 

Our article is organised as follows. In Sec. || we 
present the quantum kinetic equation and review prop- 
erties of the source term. The collision term is described 
in Sec. III. This completes the specification of the model 
and so our results appear in Sec. IV. Section ^ presents 
some concluding remarks. 



II. DISTRIBUTION FUNCTION 
A. Kinetic Equation and Source Term 

We assume that in its wake an URHIC leaves a high 
energy-density electric field, which occupies a large (in 
fact, unbounded) spacetime volume. This excited do- 
main decays via a Schwinger-likc mechanism, producing 
an unequilibrated plasma of highly energetic quarks and 
antiquarks. That system evolves and equilibrates, form- 
ing a component of the QGP. Hadronisation only takes 
place at a later stage, when the temperature and density 
of the equilibrated system fall below some critical values, 
and herein we do not consider that process. 

We represent the excited domain by a spatially homo- 
geneous, time-dependent Abelian vector potential, A^(t), 
and work in the temporal gauge: Aq — 0. The spa- 
tial part of the vector potential defines the z-direction; 
i.e., A{t) = {0,0, A{t)), and generates an electric field 
E{t) = —dA{t)/dt. This provides the input for the quan- 
tum kinetic equations whose solution describes the evolu- 
tion of the single parton distribution functions that char- 
acterise the produced partons: 



where — 
dt 



S±(p,t) + C±{p,t) , (1) 



e is the electric charge and +/— denotes bosons/fermions; 
and f± (p, t) gives the ensemble fraction of particles with 
a given momentum, p, at time t. The analogous equation 
for antiparticles, which is obtained via charge conjugation 
and must be solved simultaneously, yields f± {p, t) , the 
single antiparticle distribution function. The feedback 
generated by the motion of the partons is incorporated by 
coupling Maxwell's equa tion to Eq. (|l]) and its analogue, 
as we discuss in Sec. [IB. That also introduces a coupling 
between the equations for /, /. 

In Eq. (P, C±{p,t) is the collision term, which also 
couples the equations for / and /, and which we dis- 



cuss in Sec. [II; and S±{p,t) is the particle-antiparticle- 
producing source term: 

I /■* 

S±{p,t) - -W±{p,t,t) dt'W±{p,t,t')[l±2f±{p,t')] 

^ Jto 

(3) 



X cos 



2 dTu;{p,t,T) 



The effect of quantum statistics on the particle produc- 
tion rate is evident in the "±2 /±" in Eq. (||) (neglecting 
this term defines the low density limit) and in the differ- 
ent transition (or tunnelling) amplitudes 



W±{p,t,t')=eE{t') 



P{t,t') 
uj^pXt') \p{t,t') 



9±-l 



(4) 



where: g± — 2s± ^ 1, with s+ ~ 0, s_ = 1/2; the three- 
vector momentum p — (p*i,p||); the transverse mass- 
squared e\ = +p\; uj^{p,t,t') = e\ +p'^{t,t'); and 

p{t, t') =p\\-e [A(t) - A{t')] ^pn+e J dr E{t) . 

(5) 

(Equation (||), which describes the action of the field on 
the particles, is just a re-expression of the Lorentz force 
law: dp{t,t')/dt — eE{t).) The source term is nonlo- 
cal in time and that can be important on short time- 
scales in strong fields: if the fields are strong enough the 
time duration of a tunnelling event and the time between 
successive events, which is set by the particles' Comp- 
ton wavelengths, are similar, and the processes interfere, 
with observable consequences in the distribution func- 
tion. In addition, strong fields enhance the differences 
between fermion and boson production. These features 
were highlighted in Refs. 0, 

With our initial setup we have a simple, algebraic sym- 
metry between the particle and antiparticle distribution 
functions: 



f{p,t) = f{-p,t). 



(6) 



and hence it is only necessary to consider Eq. (^ explic- 
itly. To simplify this kinetic equation we follow Ref. 
and introduce two real auxiliary functions 



u±{p,t) 



v± (p, t) 



to 

X sin 
t 



dt'W±{p,t,t')[l±2f±ip,t')] 



dTLuip, t, t) 



I dt'W±{p,t,t')[l±2f±{p,t') 

2 j^^ dTUj{p,t,T) 



(7) 



(8) 



such that for C±(p,t) =0: u'^^ + v^T (1±2/)^ = const., 
with the initial conditions /(io) = ^^(^o) = ^^(^0) = 0. 



FIG. 1: Momentum dependence of the single parton distri- 
bution function: f( p, t = 10 fm), obtained by solving Eqs. (H)- 
( pT[ ) and (p^)-(p^, with the impulse profile in Eq. (jlj) 
{Ao = IGeV^, 176 = IGeV, and mb = 1/5; i.e., Set 4 in 
Table ^ and neglecting collisions, "transverse momentum" 
represents \p±\ and "longitudinal momentum" represents . 
The irregular structure makes clear that this is not the dis- 
tribution function of a system in equilibrium. 



This permits us to rewrite Eq. as a system of coupled, 
first-order differential equations 



(9) 

2u;{p)v±, (10) 
W±ip,t)[l±2f±]-2uj{p}u±, (11) 



where W±{p, t) and uj{p) denote, respectively, W±{p, t, t) 
and Lij{p,t,t). While this complex is qualitatively identi- 
cal to Eq. (nl), it is simpler to treat numerically. 

The kinetic equation describes pair creation for both 
bosons and fermions. In the following we identify the 
fermionic degrees of freedom with the quarks and an- 
tiquarks produced in an URHIC and restrict ourselves 
solely to the fermionic case, henceforth suppressing the 
± subscript. 



B. Maxwell's Equation and Internal Currents 

The kinetic equation depends nonlinearly on the time- 
dependent electric field. The URHIC provides the impe- 
tus for this field: it provides an external field, which we 
model via 



(12) 



This profile "switches-on" ai t ^ —2h and off at t ~ 2&, 
and attains its maximum magnitude of aX t — 0. 



FIG. 2: Time evolution of the single parton distribution 
function: f(p± = Q, p\\,t) in the low density limit. (Model 
parameters are as in Fig. hi) The regular (plasma) oscillation 
is driven by the back-reaction phenomenon. 



The external field, E^x, polarises the vacuum, generat- 
ing a polarisation current that depends on the dielec- 
tric properties of the medium, which are encoded in the 
source term, and promotes the spontaneous production of 
particle-antiparticle pairs, which it then accelerates, gen- 
erating a conduction current that depends on the particle 
distribution function. A consequence of the URHIC then 
is the appearance of a two component internal current 
and an attendant internal electric field: 

-Einit) = jm = jcond{t) + jpol{t) , (13) 

where the fully renormalised currents are p4| 



jcondit) = 2Npel ^ fip, t) , 



jpoiit) = NpC 



(27r)3 uj{p} 
d^p e± 



(14) 



(27r)3 uj{p} 



v{p,t) 



eE{t)£j_ 
4lu^{p) 



(15) 



with Np = 18, as described in Table |[ 

In the absence of collisions, Eqs. (|9|)-(pl]) and (p^)-(p^ 
with 



E{t) = EeAt) + E,n{t) ., 



(16) 



form a closed system of coupled equations whose solu- 
tion provides the time-dependent electric field and single 
parton distribution function. We can use this system 
to illustrate some of the effects we have mentioned. In 
Fig. we depict the momentum dependence of the distri- 
bution function. The irregular structure is produced by 
interference effects in the non-Markovian source term, 
which arise because the tunnelling time is of the same 
magnitude as the Compton wavelength of the produced 
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TABLE I: Parameter sets used to specify our model of an 
URHIC, Eq. (p^. They yield equilibrium values of thermo- 
dynamic quantities (three rightmost columns) that are con- 
sistent with those expected in a QGP, with Sets 3, 4 approx- 
imating RHIC-like conditions. We use mb = 1/5, a strong 
coupling: e = 1, Tc = 1, and consider 3 quark flavours, which 
explains the factor of A'j, = 18 = 2apin 3flavour 3coiour that ap- 
pears frequently. 





Ao [GeV^] 


b [fm] 


£ [GeV/fm^] 


T [GeV] 


n [fm-3] 


Set 1 


0.25 


0.2 


1.0 


0.20 


1.4 


Set 2 


0.40 


0.2 


2.5 


0.26 


3.0 


Set 3 


0.75 


0.2 


12 


0.38 


10 


Set 4 


1.0 


0.2 


24 


0.44 
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particles, and by the feedback mechanism. This struc- 
ture is averaged out in the ideal Markov limit |l^ and 
makes very clear that this is not the distribution func- 
tion of a system in equilibrium. Figure |2| depicts the time 
evolution of the calculated distribution function and the 
feedback that characterises the behaviour of the internal 
currents is manifest in the obvious, regular plasma oscil- 
lation. For comparison, in Fig. ^ we depict the distribu- 
tion function obtained after the inclusion of dissipative 




effects, to be discussed in Sec. ID 



FIG. 3: f{p, t — 10 fm) obtained with the inclusion of parton- 
parton collisions, as described in Sec. |ll4 Collisions pro- 
mote equilibration of the system, which is evident in the now 
smooth distribution function cf. Fig. 0. (Model parameters 
are as in Fig. |l|.) 



antiparticle (pa) collisions and reflect the symmetries of 
our initial conditions. Hence we employ 



III. EQUILIBRATING COLLISIONS 



Once the particles are produced they are accelerated 



by the electric field and, as we saw in Sec. II B, if eE ~ 
then large amplitude, high frequency plasma oscillations 
appear. This collective effect, which is a hallmark of our 
approach, may have observable consequences in experi- 
ments aimed at producing an equilibrated QGP. Whether 
that is the case or not can only be determined once the 
effect of parton-parton collisions is incorporated. 

The general nature of the dissipative collision term is 
known: it too is non-Markovian and can produce parti- 
cles 
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However, its complexity mitigates against its 
use in semi-quantitative, exploratory studies and hence 
herein we employ a simple RTA |l|, |l|, |l|, |l9[ |o[ |2^, 

MM-, 

In this approach the detailed description of parton- 
parton scattering is replaced by a continuous viscosity 
term, which involves a time-dependent parameter that 
is identified with the collision period. In addition we 
suppose that the thermodynamic laws are valid at each 
time t, and this assumption of local-equilibrium provides 
for an internally consistent definition and calculation of 
time-dependent thermodynamic variables, such as tem- 
perature and energy density. Of course, accepting a phys- 
ical interpretation of these quantities only makes sense 
once the violent effects of the URHIC have subsided and 
the quantities are evolving slowly with t. 

The collision term in the kinetic equation for 
f{p, t) must describe particle-particle {pp) and particle- 



C{p,t;T,u-)^ 

^ [f,,{p,t;T{t),um))- f{p,t)] 



Tpp{t) 
1 



Tpa{t) 



[f,,{-p,t;Tit),u'^it))~f{p,t)] (17) 



where Tpp{t) — Tpa{t) — T(t) is the time-dependent relax- 
ation time and 



feq ip,T, u") 



exp 



T{t) 



(18) 



is the quasi-equilibrium distribution function, with p^ the 
quarks' 4-momcntum. The other quantities in Eq. (pT 
are: 



U^it) = (l,0,0,u(t))[l-M(t)2]-(l/2) 



(19) 



the hydrodynamical four- velocity; and T{t), the local- 
equilibrium temperature, both of which we define below. 

Our RTA encodes all the complicated effects of parton- 
parton collisions in a single time-dependent quantity, the 
relaxation time: T{t). It is a measure of the time between 
successive collisions and as such can be expressed as |23] 



Tit) = T, ^ , 
v{t) 



(20) 



where \{t) = \/[n{t)Y/'^ is the mean interparticle sepa- 
ration, since 



<t)-Np i^fm 



(21) 
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is the mean particle number density, and 
v{t)^\p\f{t)/ef{t) 



(22) 



is an average speed. (See Eqs. (g3|) and (|27|).) The only 
parameter in our implementation of the relaxation time 
approximation is then Tc, the dimensionless constant of 
proportionality. 

Returning to Eq. (^7|) , we define the local-temperature 
by requiring that at each time, the mean particle en- 
ergy density in the plasma is identical to that in an equi- 
librated plasma at a temperature T{t); i.e.. 



^u;{p)[f{p,t)-Zj{p,t)] (23) 



(2^ 



(24) 



Similarly, u{t) in Eq. ([l9| ) is defined via the requirement 



Pf{t) = N, 



(Pp 
^ / (2^ 

d?p 



p[fip,t)~zf{p,t)] (25) 
Pfeg{p,t); (26) 



i.e., that the mean particle three-momentum is the same 
as that of an equilibrated plasma characterised by an 
hydrodynamical particle velocity u(t). The average mag- 
nitude of the momentum is 



\p\f{t)^Np 



(Pp 
(2^ 



\j^[f{p,t)-zf{p,t)]. (27) 



The new element in these equations, 



Zfip.t) = 



e £_L 
4w3(p) 



2 r 



(28) 



is a regularising counterterm, determined via the same 
proc edur e jl^ that yields the renormalised currents in 
Sec. [IB , which ensures that the integrals involving the 
calculated distribution function are finite. This countert- 
erm itself exhibits "memory effects;" i.e., it is sensitive 
to the time-history of the electric field. 

We judge that using a time-dependent relaxation time 
is an improvement over previous work that used T(i) — 
const.; e.g., Refs. |l5|, because it is manageable and 
better models the conditions produced by an URHIC. 
Clearly, just after the impact the parton number density 
is small and hence the time between successive collisions 
is large. Over time, however, the density of produced par- 
tons increases, leading to a reduction in the interval be- 
tween collisions. These features are crudely reflected by 
the evolution of the relaxation time described in Eq. (20). 
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-1 1 Time [fm] 345 

FIG. 4: Time evolution of the total electric field. The effects 
of feedback are evident in the oscillatory behaviour of E{t). 
The viscous collision term damps these oscillations in a char- 
acteristic time Tc/m ^ 1 fm. The external impulse electric 
field is also depicted. (Parameters: Set 2, Table |l[) 




Time, [fm] 



FIG. 5: Time evolution of the calculated energy density, 
Eq. ( p3[ ) . It is negative in the neighbourhood oi t — because 
of the large vacuum polarisation produced by the strong ex- 
ternal field (the URHIC). This feature is quickly compensated 
by particle production in the aftermath of the impact. Col- 
lisions subsequently yield an equifibrated QGP. (Parameters: 
Set 2, Table |.) 



The final form of our collision term is thus obtained 
from the combination of Eqs. ( p7| ) and (|20|); i.e., 

C(p i) = M ( feMT{t),u{t)) - fip,t) 



feqi-p,Tit),uit))-f{p,t) 



(29) 



It provides additional nonlocal feedback in the solution 
for f{p,t). We emphasise again that our RTA is based 
on the assumption of local-equilibrium, which is valid for 
\t\ 3> b, where b is the time duration of the URHIC. For 
\t\ ^ b, however, it is of questionable validity and may 
lead to model-dependent artefacts, the misinterpretation 
of which one must guard against. 
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FIG. 6: Time evolution of the particle number density, 
Eq. (0). The equihbrium value: 3.0 fm^'^, is that listed in 
TableH (Parameters: Set 2, Table |.) 
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FIG. 7: Time evolution of the entropy/particle, Eq. (pO[), 
with and without collisions. It is intuitively obvious why the 
entropy is larger when particle-particle collisions are incor- 
porated: in equilibrating the system, collisions increase the 
degree of disorder. (Parameters: Set 2, Table ffl) 



IV. NUMERICAL RESULTS 
A. Thermodynamic Parameters 

Our results are obtained via the simultaneous solution 
of Eqs. (Ol)-®' (il)-®, using the colhsion 

term in Eq. (|29|), which we accomplish using a fourth- 
order Runge-Kutta procedure. The solution, /(p, t), fully 
describes the plasma's evolution, from its creation to 
equilibration. In addition we obtain the time-dependent 
vector potential, electric field and currents, and also the 
quasi-equilibrium temperature and collective velocity. 

To explore the solution's properties we have employed 
a range of parameter sets, which are listed in Table |, 
and in Figs. we demonstrate the effect of collisions 
by comparing the solution obtained using Tc ~ 1.0 with 
that obtained in the coUisionless limit Tc oo. (NB. As 
evident in the figures, the plasma period satisfies rriTpi w 
2.5 ~ Tc- For Tc <tC TfiTpi plasma oscillations are not 
observable 113]; i.e., the system is overdamped.) 



FIG. 8: Time evolution of the quasi-equilibrium tempera- 
ture. The equilibrium value reached here: 0.26 GeV, is greater 
than the anticipated critical temperature for QGP formation. 
(Parameters: Set 2, Table ^) 



In Fig. H we see that the URHIC generates a strong 
electric field, which produces particles that sustain the 
field for a time that depends on the collision frequency: 
a large value of Vc = 1/tc means a short-lived electric 
field and rapid equilibration. 




Time [fm] 



FIG. 9: Time evolution of the collective velocity, u(€). 
It's behaviour clearly signals the collective plasma oscillation 
that is characteristic of the flux-tube production mechanism. 
The plasma oscillation is damped by collisions. (Parameters: 
Set 2, Table |.) 

Figure [q d epicts the time evolution of the energy den- 
sity, Eq. (|23|). It is negative on a small domain around 
t = because of the large vacuum polarisation induced 
by the strong external field. This instability is quickly 
corrected by rapid particle creation' a correlation that is 
apparent in a comparison of Fig. ^ with Fig. which 
portrays the particle number density, Eq. (|2l|). The den- 
sity reaches a higher value in the absence of collisions 
because the field-current feedback allows unhindered, re- 
peated bursts of particle-pair creation. The Pauli prin- 
ciple does not significantly retard the process because, 
while the particles are preferentially produced with small 
momenta, the field rapidly accelerates them. As Table | 
shows, the number density increases with increasing Aq\ 



FIG. 10: Time evolution of the dilepton production rate. 
The plasma oscillation is evident in this collisionless case, as 
is the feature that the dileptons are preferentially produced 
at low-M^. (Parameters: Set 4, Table |.) 



FIG. 11: Time evolution of the dilepton production rate. 
Frequent collisions rapidly damp the plasma oscillation, forc- 
ing the production rate to settle at a constant value. (Param- 
eters: Set 4, Table ^) 



i.e., with increasing impulse field strength. 

In Fig. 1^ we plot the behaviour of the entropy /particle, 
where the entropy density is: 



sit) = -Np 



(2^ 



/(p,t)ln/(p,t) 



(30) 



The calculated quasi-equilibrium temperature is de- 
picted in Fig. ||. After rising quickly, it settles into a 
slow evolution once the external field, Eq. (p^), has sub- 
sided. This marks the beginning of the domain on which 
the concept of local equilibrium is valid. 

The local velocity is plotted in Fig. ||. It shows, as 
one would intuitively expect, that the produced particles 
are accelerated, reaching their maximum velocity when 
the electric field vanishes (cf. Fig. then decelerated as 
the field reverses direction. They stop, and then reverse 
direction and are accelerated by the reversed field to a 
new maximum velocity. The repetition of this pattern is 
the collective plasma oscillation, which is clearest in this 
figure. Of course, in equilibrating the system, collisions 
act to destroy the pattern. 



B. Dilepton Production 

The plasma oscillation is evident in the electric field 
and in each of the thermodynamic variables but none of 
these quantities are directly measurable. Is there any 
way that this characteristic signature of the flux-tube 
production mechanism can be observed? 

The thermal dilepton production rate may provide a 
means. Dileptons produced in the URHIC do not in- 
teract strongly and hence those produced soon after the 
impact carry and transmit information about the pre- 
equilibrium stage of the plasma. The dilepton production 
rate from our (quasi-)equilibrium three quark system can 



be estimated using | 
dN 



37r3 



(l-4mf/Ar)5 1 + 2 



AP 



M4 



poo 

X / deide2f{ei)f{e2)e{AP-Ml)9{Ml-A'P), 

J m 

(31) 

where: a — 1/137; m; is the lepton mass, m the quark 
mass and we use mi = m; is the invariant mass of 
the produced dilepton pair; 



All = 2m2 + 2eie2±2|pi||p2|; 



(32) 



and the time-dependence of the distribution fimctions is 
implicit. 

We obtain the distributi on f unctions, /(ei,2) in 
Eq. (|3l|), as described in Sec. IV A , using the RHlC-like 



parameter-Set-4 in Table |. Knowing them, the calcula- 
tion of the dilepton production rate is straightforward. In 
Fig. |l^ we display that rate as a function of (A'P , t) . The 
plasma oscillation is evident in the time evolution, with 
more dilepton pairs being produced when the electric 
field is strongest and the parton production rate peaks, 
and no pairs being produced when the electric field van- 
ishes. The effect of collisions is to drive the system to 
equilibrium where the thermal dilepton production rate 
becomes constant, as shown in Fig. |ll|. 

Clearly, the plasma oscillations generate a signal. How- 
ever, the time evolution of the dilepton production rate 
is a difficult quantity to measure. An easier quantity is 
the time-integrated rate: 



Pi+i-{t) ■■= 



dN 



d^xdAP 



= / dt' 



dN 



dt' d^x dAP 



(33) 



Does a signal survive in this observable? 
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FIG. 12: Time integrated dilepton production rate, Eq. (p3[), 
calculated without collisions and hence in the presence of a 
persistent plasma oscillation. (Parameters: Set 4, Table 0.) 




FIG. 13: Time integrated dilepton production rate, Eq. (^sj), 
calculated under the influence of frequent collisions that 
rapidly equilibrate the plasma. (Parameters: Set 4, Table ^) 

In Fig. |l^ we plot Pi+i-{t) obtained in the absence 
of collisions, which must be compared with the function 
in Fig. |l^ that was obtained with the inclusion of colli- 
sions. Our simple model yields rates that are comparable 
with other estimates; e.g., Ref. |2^, and the comparison 
of the figures shows that plasma oscillations generate an 
enhancement in the number of dileptons, Pi+i-(t), which 
is as large as a factor of 2 at t = 5 fm, given the Set 4 
initial conditions in Table |. (Note that collisions elimi- 
nate the plasma oscillation so Fig. |l^ can be thought of 
as the oscillation- free scenario.) 

A further illustration of the effect is provided in Fig. |lj. 
The equilibrium energy per particle, s/n, grows with the 
violence of the collision; i.e., with the value of Aq, as do 
the amplitude and frequency of the plasma oscillations. 
This effect is responsible for the sharp increase in the 
dilepton yield, evident in Fig. |l^, for the most energetic 
collision in Table ^. At lower values of e/n the effect 
of plasma oscillations is suppressed by collisions, but for 
the Set 4 initial conditions the magnitude and frequency 



of the plasma oscillation are large enough to make their 
action evident in spite of the damping, at least in our ide- 
alised treatment. (NB. A rapid expansion of the plasma 
will alter the initial conditions required for plasma oscil- 
lations to have observable consequences.) 



s.o 



E 4.0 

II 

is 



0.0 
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, I — . — . — . — , — . — . — . — , — . — . — . — , — . — . — . — ^ 
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FIG. 14: Pi+i- {t = 5fm) at the low invariant mass, , for 
which it takes its maximum value, as a function of the equi- 
librium energy /particle attained in the collisions described 
by the parameter sets in Table |. The sharp increase in Pi+i- 
occurs after the plasma oscillation overpowers the viscous col- 
lision term. 



V. EPILOGUE 

We have used a quantum kinetic equation coupled with 
Maxwell's equation to explore the formation and equili- 
bration of a strong field plasma. Using a simple impulse 
model for the URHIC, which produces RHIC-like con- 
ditions, we find that the non-Markovian aspects of the 
source term do not generate observable effects. However, 
the field-current feedback, which is characteristic of the 
production of strongly coupled charges by a strong field, 
manifests itself in the appearance of plasma oscillations 
in the thermodynamic observables. The oscillations are 
also evident in the production rate of thermal dileptons 
and, while the time evolution of this rate may not be 
measurable, the plasma oscillations act to significantly 
enhance the time integrated rate. The effect is marked 
by a sharp increase in the dilepton yield when the en- 
ergy per particle becomes large enough to generate a high 
frequency and large amplitude plasma oscillation, which 
initially overwhelms the effect of collisions. 

While the magnitudes of the quantities we calculate 
are phenomenologically reasonable, the primary results 
of our study are qualitative and many improvements 
are possible. Our relaxation time approximation to the 
collision term is an intuitive and practical tool but a 
more realistic connection with the actual collision process 
would provide a systematic and well-constrained quanti- 
tative improvement. A simpler step is the introduction 
of a strongly momentum-dependent dressed-parton mass, 
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which is an essential feature of QCD |28|. That can lave 
a significant impact on the evolution of the plasma, pro- 
moting plasma oscillations |2^, and also on its si bse- 
quent hadronisation ||30| . Perhaps the most significant 
defect of our study is the use of an Abelian model fo ■ the 
colour fields and progress with a non- Abelian transport 
equation would be a marked improvement [n9| . 



ACKNOWLEDGMENTS 



Ti mdy 



We thank R. Alkofer, V. Morozov and P.C 
for helpful discussions. This work was supportec 
the Deutsche Forschungsgemeinschaft under project 
SCHM 1342/3-1 and 436 RUS 17/102/00; the US Deiart 
ment of Energy, Nuclear Physics Division, under contract 
no. W-31-109-ENG-38; the US National Science 
dation under grant no. INT-9603385; and partly bj 
Russian Federation's State Committee for Higher 
cation under grant no. N 97-0-6.1-4. This research 
fited from the resources of the National Energy Research 
Scientific Computing Center. 



by: 
nos. 



F Dun- 



le: 



REFERENCES 

[1] E. Laermann, Fiz. Elem. Chastits At. Yadra 30 720 
(1999) [Phys. Part. Nucl. 30, 304 (1999)]. 

[2] CD. Roberts, and S.M. Schmidt, Prog. Part. I|Jucl. 
Phys. 45, SI (2000). 

[3] S.A. Bass, M. Gyulassy, H. Stocker and W. Gre 
J. Phys. G G 25, Rl (1999); S. Scherer et ai, I 
Part. Nucl. Phys. 42, 279 (1999); U. Heinz 
M. Jacob, "Evidence for a new state of matter 
assessment of the results from the CERN lead h 



Alm- 
ner, 



programme," |nucl-th/0002042| . 
[4] L.V. Gribov and M.G. Ryskin, Phys. Rept. 18f , 29 

(1990) . 

[5] B. Andersson, G. Gustafson, G. Ingelman antl T. 

Sjostrand, Phys. Rept. 97, 33 (1983). 
[6] X. Wang and M. Gyulassy, Phys. Rev. D 44, ^501 

(1991) . 

[7] K. Geiger, Phys. Rept. 258, 237 (1995). 
[8] B. Andersson, G. Gustafson and B. Nilsson- 

quist, Nucl. Phys. B 281, 289 (1987); K. We: 

Phys. Rept. 232, 87 (1993). 
[9] M. Gyulassy and T. Matsui, Phys. Rev. D 29, 419 

(1984); K. Kajantie, R. Raitio and P.V. Rm.ska 

nen, Nucl. Phys. B 222, 152 (1983); G. Baym, 

B.L. Friman, J. P. Blaizot, M. Soyeur and W. Czyi, 

Nucl. Phys. A 407, 541 (1983). 
[10] Y. Kluger, E. Mottola, and J.M. Eisenberg, Ifhys. 

Rev. D 58, 125015 (1998). 
[11] S.M Schmidt D. Blaschke, G. Ropke, S.A. SmoJyi 

sky, A.V. Prozorkevich and V.D. Toneev, In 

Mod. Phys. E 7, 709 (1998). 
[12] S.A. Smolyansky, G. Ropke, S.M. Schmidt 

Blaschke, V.D. Toneev and A.V. Prozorkevich, |hep- 



thc 
du- 
ne- 



rog. 
and 
An 



an- 
J. 

D. 



[13] 

[14] 

[15] 
[16] 

[17] 

[18] 

[19] 

[20] 

[21] 
[22] 
[23] 

[24] 



Bloch 



r iS, 



ttc 



[25] 

[26] 
[27] 
[28] 

[29] 
[30] 



ph/9712377 
S.A. Smoly 
pairs in 
Nonequilibr 
27-30, 1998 
S.M. Schmi 
zorkevich, £ 
Rev. D 59 
J.C.R 
CD. Robe 
Vinnik, Phj 
J.C.R Blocl 
Rev. D 61 
Y. Kluger, 
and E. Mo 

F. Cooper, 
and B. Svet 
K. Kajantie 
(1985); G. 
Phys. Rev. 
B. Banerj( 
Phys. Lett 

G. C Nayak 
6877 (1997) 
A. Bialas 
Nucl. Phys 
S. Mrowczh 
G. Baym 
G.C. Nayak 
Greiner, 
at RHIC 
L.D 

545 (1985); 
and S. Prem 
A. Drees, 
edited by D 
(World 
T.S. Biro 
(1997); W 
Eur. Phys. 
Wilk and Z 
CY. Wong, 
Collisions, 
D. Dutta, 
Choudhury, 
CD 

ment and 



S.M. Schmidt, A.V. Prozorkevich, and 
msky. Creation of boson and fermion 
strong fields, Proceedings 'V. Workshop on 
urn Physics at Short Time Scales', April 



W. 



Phys 



and 
McLerran 



Scieitific, 
aid 



M. 



Roberts 



nucl-th/000 ro54 



hep-ph/9809233 



D. Blaschke, G. Ropke, A.V. Pro- 
A. Smolyansky and V.D. Toneev, Phys. 
094005 (1999). 

V.A. Mizerny, A.V. Prozorkevich, 
S.M. Schmidt, S.A. Smolyansky, D.V. 
s. Rev. D 60, 116011 (1999). 

CD. Roberts and S.M. Schmidt, Phys. 
117502 (2000). 
M. Eisenberg, B. Svetitsky, F. Cooper 
ola, Phys. Rev. Lett. 67, 2427 (1991); 
J.M. Eisenberg, Y. Kluger, E. Mottola 
itsky, Phys. Rev. D 48, 190 (1993). 
and T. Matsui, Phys. Lett. B 164, 373 
Gatoff, A.K. Kerman, and T. Matsui, 
3 36, 114 (1987). 

R.S. Bahlerao and V. Ravishankar, 
B 224, 16 (1989). 

and V. Ravishankar, Phys. Rev. D 55, 
ibidC 58, 356 (1998). 

Czyz, A. Dyrek and W. Florkowski, 
B296, 611 (1988). 
3ki, Phys. Rev. D 39, 1940 (1989). 
s. Lett. B 138, 18 (1984). 
A. Dumitru, L. McLerran and W. 
Ejquilibration of the gluon-minijet plasma 
LHC," |hep-ph/000l"20^ . 

and T. Toimela, Phys. Rev. D 31, 
D. Boyanovsky, H.J. de Vega, R. Holman 
Kumar, Phys. Rev. D 56, 5233 (1997); 
Jnderstanding Deconfinement in QCD, 
Blaschke, F. Karsch and CD. Roberts 
Singapore, 2000) pp. 285-290. 
C. Greiner, Phys. Rev. Lett. 79, 3138 
Alberico, A. Lavagno and P. Quarati, 
J. C 12, 499 (2000); O.V. Utyuzh, G. 
Wlodarczyk, J. Phys. G 26, L39 (2000). 
Introduction to High-Energy Heavy-Ion 
(|World Scientific, Singapore, 1994). 
K. Kumar, A.K. Mohanty and R.K. 
Phys. Rev. C 60, 014905 (1999). 

"Continuum strong QCD: Confine- 
lynam ical chiral symmetry breaking," 
P. Maris, "Continuum QCD and 



light mesonj^," |nucl-th/0009064 . 
A.V. Prozc 
Roberts anc 
ma Oscillat - 
ter in Astrf 
many, Nov 

P. Bozek, Yl.B. He and J. Hiifner, Phys. Rev. C 57, 
3263 (1998), 



zo]tkevitch, D.V. Vinnik, M.B. Hecht, CD. 
S.M. Schmidt, Pair Creation and Plas- 
ons, in the Proceedings of "Quark Mat- 
and Part iclephysics," Ros tock, Ger- 



7-29, 2000, |aucl-th/0012039 



